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Abstract 


The gelation kinetics of silica nanoparticles is a central process in physical chemistry, yet 
not fully understood. Gelation times are measured to increase by over four orders of mag¬ 
nitude, simply changing the monovalent salt species from CsCl to LiCl. This striking effect 
has no microscopic explanation within current paradigms. The trend is consistent with the 
Hofmeister series, pointing to short-ranged solvation effects not included in the standard col¬ 
loidal (DLVO) interaction potential. By implementing a simple form for short-range repulsion 
within a model that relates the gelation time-scale to the colloidal interaction forces, we are 
able to explain the many orders of magnitude difference in the gelation times at fixed salt con¬ 
centration. The model allows to estimate the magnitude of the non-DLVO hydration forces, 
which dominate the interparticle interactions at the length-scale of the hydrated ion diame¬ 
ter. This opens the possibility of finely tuning the gelation time-scale of nanoparticles by just 
adjusting the background electrolyte species. 

Keywords: gelation time; Hofmeister series; silica nanoparticle gel; hydration forces 


The total interaction potential in acqueous suspensions of charged colloidal particles is often 
taken as the sum of the van der Waals attraction and a simple approximation of the electrostatic 
double-layer repulsion, forming the classical Derjaguin-Landau-Verwey-Overbeek (DLVO) poten¬ 
tial^. The electrostatic repulsion can be overcome by adding salt, thus increasing the screening of 
electrostatic repulsion, and lowering the energy barrier against aggregation^^. Colloid particles then 
typically aggregate into clusters, which grow over time with the possibility of forming a sample- 
spanning network. Silica particles, and specifically Ludox, have been a classical model system 
for colloidal physical chemistry, and also have a key place in industrial processing, coatings, ink 
receptive papers, metal casting, refractory products, and catalysts. 

DLVO theory predicts that there is a salt concentration at which the suspension aggregates 
(critical coagulation concentration); here the maximum DLVO interaction and its derivative are 
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Figure 1: The time required for gelation of the sample is related to the interparticle interaction 
potential, and the specific hydration of the ions has a huge effect, controlling the inter-colloid 
hydration repulsion at a /z < A. (a) Schematic, and (b) plot, of the “total” interaction obtained 
in this work, including the strongly ion-specific short-ranged repulsive hydration potential. Also 
illustrated are the classical terms in the DVLO interaction. At very short distances, on the order of 
the hydrated ionic diameter, the repulsive shoulder dominates the interaction, (c) The macroscopic 
sample initially flows as the vial is inverted (left), while after some time it becomes solid (right). 
This sample has <)>Ludox = 0- 140 and 374 mM NaCl. 


both zero (i.e. there is no energy barrier against aggregation). This concentration is proportional 
to where z is the valency of the salt ions, and is one of the successes of DLVO theory!^. It is 
well known that DLVO theory breaks down completely for high salt concentrations (above O.IM, 
which unfortunately is the regime of biological interest)®^. In this range, the basic assumptions 
of point charges, a solvent continuum and neglecting ion-surface adsorption and dispersion forces, 
are called into question. On a fundamental level, even the assumption that electrostatic forces and 
dispersion forces are additive is incorrect!^. 

It is shown in this letter that there is a spectacular failure of DLVO theory in estimating gelation 
times, for identical particles, in the presence of different monovalent salts even at low concentra- 
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tion. Changing the salt type dramatieally affeets the aggregation proeess. Using a model for 
relating the gelation kineties to partiele interaetions, a non-DLVO hydration repulsion is eharae- 
terised, and the dramatie ehanges in gelation kineties are explained mieroseopieally in terms of ion 
solvation and its interplay with the eharged eolloid surfaee. The proposed framework will make 
it possible to finely tune the gelation rate of nanopartieles simply by the ehoiee of monovalent 
eleetrolyte speeies in the eolloidal solution. 

A minimum of eontext and eoneepts proposed to explain the Hofmeister series are useful to 
the reader. It was first shown by Hofmeistei® that the stability of a eolloidal solution (he made 
observations on proteins, whieh were then investigated by others^^, while other work investigated 
eolloidal partielesl^l^^^ ean be drastieally different upon the addition of different salts of the same 
valeney, even if all other parameters (sueh as the salt eoneentration) are kept eonstant. Eleetrolytes 
eould be arranged aeeording to their effieieney in salting out protein (now known as “Hofmeister 
series”). The effeet is understood to be related to how salt ions strueture the water around them¬ 
selves. For monovalent eations the series is NHJ, Cs'*', Rb'*', K'*', Na'*', LU, from most ehaotropie 
(weakly hydrated, strueture breakers) to most eosmotropie (strongly hydrated, strueture makers), 
and whieh extreme is most destabilising depends on the surfaee properties of the eolloids^. 

There are numerous, partially eonflieting, theories as to the origin of these short-ranged ion- 
speeifie eolloidal interaetions, linking to the ioni sizd ^^b4|i 5| Strongly polarisable ions are large 
and have more diffuse eleetron elouds. The energy penalty for being less well hydrated (for exam¬ 
ple due to adsorption at an interfaee) is low for sueh ions beeause the eharge ean be easily redis¬ 
tributed^. The deereasing size trend in going from ehaotropie to eosmotropie in the Hofmeister 
series is eonsistent with this pietureff^. However the effeetive polarisability of the ions eonsists of 
eontributions from both the ion itself and the solvent moleeules in its hydration shell^. This by 
extension may signifieantly augment the dispersion forees, and give further ion-speeifie interae- 
tions^^. Another suggestion is that the large eleetrie field on the eolloidal surfaee, arising from the 
finite size of the eounterions, results in the ions aequiring appreeiably large effeetive polarisabili- 
tiesHH. 
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Figure 2: There is a very strong power law dependence (approximately -6 exponent) between gela¬ 
tion time and salt concentration, and a striking difference between the five monovalent salt species. 
The experimental observations (markers) are well recapitulated by the theoretical predictions (one 
parameter fits, as described in the text) (solid lines). Data is obtained from samples with volume 
fractions 0.13 and 0.14, and are undistinguishable. For each curve, the only free fit parameter is the 
hydration force amplitude Fq, which is a function dependent on the salt-concentration (see Fig.[^. 
The theory is calculated assuming volume fraction 0.133, and provides a match with the data for 
values of Fq well within the typical range of 10^ to 5 x 10^ Nm^^ (retfi^. 

Experiments with (negatively charged) mica surfaces showed that there is adsorption of cations. 
More hydrated cations (such as Li'*') are adsorbed only at high salt concentrations, while the less 
hydrated ions adsorb at lower concentrations. However, once adsorbed, the cosmotropic ions retain 
part of their hydration layer. This gives rise to a repulsive interaction as two surfaces approach 
each othei^. For mica, for example, cosmotropic ions are thus much more efficient at providing 
stabilisation than chaotropic ions. 

Models have been made to describe ion-specific distribution of ions near surfaces, and their 
surface adsorption: an important factor is the ion diameter, and whether or not the ions are hy¬ 
drated^. Chaotropic ions have smaller effective diameters (as they are not hydrated), and they 
can adsorb to the surface. Based on this model it was possible to calculate the critical coagulation 
concentration for a range of salts^. 

A fully quantitative description of hydration interactions should include all factors outlined 
above (finite ion sizes, discrete nature of the solvent, many-body dispersion forces and polarization 
effects); Such a description does not exist, and is beyond the scope of this letter, which instead aims 
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to provide motivation, guidance and useful constraints for future models. We take a simplified 
approach to modeling the hydration interactions, which is described below, and we show how 
the total interparticle interaction can be related to the gelation times, and thus measured. The 
DVLO potential Vdvlo is a linear superposition of an attractive van der Waals potential Vvdw> 
an electrostatic repulsion potential Vr and a short-ranged Born repulsion potential V^, which we 
neglect in this study due to the fact it does not impact on the gelation process. The attractive 
component is given by 



( 1 ) 


where a is the colloid radius and r = 2a + his the colloidal centre-to-centre separation, see Fig. [^a) . 


The Hamaker function A// can be written in the forml^A//(r) = Ae=o/icr(^) +^e>o//-ef(^)>where 


Ae=o is the zero frequency contribution which is screened by the counterions through the screening 
function fserif), and Ae>o is the non-zero frequency contribution which is mitigated by retardation 
through the retardation function fretir) (their full form is given in SI for completeness). 

The DLVO theory is based on a number of assumptions, the most important of which is the 
linearized Poisson-Boltzmann treatment of the electric-double layer repulsion, which is valid only 
within the Debye-Hueckel limit for the surface potential, i.e. for potentials lower than 20-25 mV. In 
our calculations, however, we used an extension due to Sader, Carnie and Chan^^, which extends 
the validity of DLVO theory to much higher potentials: 



( 2 ) 


where Co is the permittivity of free space, £„ is the relative permittivity of water, kg is the Boltz¬ 
mann constant, e is the counterion charge and T is the temperature. The function Y(r) is given in 
full in SI, and depends on the surface potential y/g. 

The other assumptions of the theory are the following: (a) The ions are treated as point-like 
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(their finite volume and exeluded-volume effects are neglected); (b) Spatial correlations among 
ions are neglected; (c) Ion-adsorption on the colloid surface is neglected; (d) Dissociation equi¬ 
libria between charged species on the colloid surface and ions in solution are neglected. Many 
studies have shown that the interparticle potential deviates significantly from DLVO theory below 
a surface-to-surface separation h of about 2 nm in watei^. An additional repulsive potential has 
been postulated to arise from the hydration of the water due to the presence of counterions and/or 
on strongly hydrophilic surfaces. This potential is still not fully understood microscopically and, 
as outlined above, there are different competing theories relating to its origins exist in the literature. 
Since the 1970s there has been a general agreement that the effective hydration potential decreases 
exponentially from the surface!^, and can thus be taken to have the following general form: 


14 = FoTTaA^exp 



(3) 


where the fitting parameters Fq and A control the magnitude and the decay of the potential, re¬ 
spectively. This expression has been used quantitatively to successfully describe the huge energy 
barrier contributed by repulsion between structured water layers on hydrophilic surfactant-coated 
colloids, a big effect which cannot be explained by DLVO-theory alone^. 

The full potential we consider is then a sum of dispersion, electrostatic and surface hydration 
terms: 

Vtot = Kdw + 34 + 34- (4) 

The interaction potential and gelation time are linked: A recent theoretical study^has estab¬ 
lished that the gelation time tgei can be evaluated according to the following expression: 

~ 2k,[(l+c)(^)o/2]+/3’ 

where df is the fractal dimension of the clusters (which for reaction limited aggregation RLCA 
is df = 2.1^^, and ^0 = (4/3);ra^no is the volume fraction of the colloids where uq is the total 
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number of colloidal particles per unit volume, uq = N/V. The parameter c equates to c = (1 — 
0c)/0c, where 0c is the critical volume fraction at which the systems gels (the zero-shear viscosity 
diverges), which for spherical-like clusters is 0c ~ 0.64. The characteristic aggregation rate kc is 
given by the relation 


kc = 



( 6 ) 


where kagg is the rate constant of aggregation. The stability ratio W is defined as^^^^^ 


1T = 




exp(/3ytot) 
{2a + h)^G{h) 


(7) 


where Vtot is the total interparticle potential, ks is the Smoluchowski diffusion limited aggregation 
rate ks = {%/'i)kBT/fi with /i the solvent viscosity. The hydrodynamics of two spheres approach¬ 
ing is given by G{h) = (6(h/a)^ -l-4(h/a))/(6(h/a)^ -|- 13(h/a) + 2). Combining Eqs. ?? and ??, 
and rewriting the colloidal concentration hq in terms of the volume fraction 0o and the volume of 
one particle Vp, we can recover the expression 


40okfir 

WpWji' 


( 8 ) 


This identity can be inserted into Eq. ?? to obtain an explicit form for the gelation time as a func¬ 
tion of the sample material characteristics and interaction parameters: 


tgel 


SVpWil 

80okBr[(l+c)0o/2]^//3' 


(9) 


Gelation times of various samples, varying salt type in Eig. salt concentration and Eudox 
concentration in Eig.[^ were determined by checking when macroscopic samples no longer flowed. 
Clearly, there is a progressive shortening of gelation times switching samples with the same salt 
concentration, going from EiCl to NaCl, KCl, RbCl and CsCl, and this is in agreement with the 
Hofmeister series. There axe, four orders of magnitude of difference in the gelation times - quite 
remarkable! - for the same concentrations of monovalent salts. 
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Figure 3: The power law dependence of gelation time with salt concentration, and the strongly salt- 
specific gelation times (shown here are (a) NaCl and (b) KCl), are seen in samples with varying 
Ludox concentration. Values of the best fit power law exponents are listed in Table 


There is a clear power law dependence of the gelation times on the salt concentration, with 
exponents (see Table all close to —6. There appears to be a decrease of power law exponent 
(to a more negative value) as the cation becomes more chaotropic. For samples with NaCl, the 
power law exponent decreases (becomes more negative) as the Ludox concentration decreases, but 
no clear trend is observed for samples with KCl. Reerink and Overbeekl^ found similar power 
law behaviour for Agl colloids, with power law exponents from around —6 to —11, depending on 
the particle size and surface potential. Other power law exponents had been reported in classical 
literaturel^, ranging from —2 to —12. The value of —6 falls into this range, and seems very robust 
in our silica colloid data. 


A power law relation is also expected if the Ludox volume fraction is varied, at fixed salt con- 
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Table 1: Power law exponents for gelation times, from the tg^\ vs. salt concentration data in 
Figs. and 1^ error quoted is the 95% confidence interval. 


Salt 

Volume fraction Ludox 

Exponent 

KCl 

0.140 

-6.27 ±0.24 

KCl 

0.105 

-6.42 ±0.34 

KCl 

0.070 

-5.92 ±1.06 

KCl 

0.035 

-7.03 ±0.42 

NaCl 

0.140 

-5.67±0.11 

NaCl 

0.105 

-6.06 ±0.30 

NaCl 

0.070 

-6.47 ±0.56 

NaCl 

0.035 

-6.60 ±0.34 

LiCl 

0.133 

-5.49 ±0.07 

NaCl 

0.133 

-5.77 ±0.10 

KCl 

0.133 

-6.68±0.15 

RbCl 

0.133 

-7.33 ±0.06 

CsCl 

0.133 

-8.44 ±0.10 


centrationl^. If the colloidal aggregation is taking place with fractal dimension df = 2.1 (as is the 
case in the RLCA regime), then the exponent in this plot is expected to be — (J//3 + 1) = —1.7; 
the data in Fig.|^show a good agreement with this theoretical expectation. Note however that since 
(as explained later) we cannot assume to have a constant hydration force, we can’t use the data of 
Fig.|^to robustly go backwards and extract the fractal dimension. 


The theoretical framework outlined above, along with the experimental data on the gelation 
times at varying salt concentration, now allows us to obtain the hydration parameters of the mono¬ 
valent salts LiCl, NaCl, KCl, RbCl and CsCl. The first question to address is how to define the 
exponential decay length A in eq. ??. Throughout the literature this parameter has been varied 
within the range 0.2 — 1.0 nm^^ for different colloidal systems. It seems reasonable to us to set A 
as the characteristic hydration diameter of the counterions ( Cs"*^ = 658 pm, Rb+ = 658 pm, K+ = 
662 pm, Na+ = 716 pm, Li+ = 764pm)l^. 


The second, more delicate question, is how to determine the two other unknown interaction 
parameters, which are the amplitude Fq of the hydration interaction (in eq. ??), and the surface 
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Ludox volume fraction 


Figure 4: Gelation times also depend on the volume fraction 0Ludox of Ludox particles: (a) KCl 
and (b) NaCl. The solid line is a guide to the eye, illustrating the slope from a power law with 
exponent -1.7, which is expected for a fractal dimension of 2.1. 


potential i/Aq (in eq. ??). 


Our first approach was to fix the hydration force constant Fq equal to a reasonable value from 
the literature^, and allow the surface potential to vary as a function of the salt concentration. In¬ 
deed, it might be expected that the association of counterions with silica surface groups will di¬ 
minish the magnitude of the surface potential with increasing salt concentration, thus reducing the 
electrostatic repulsion and speeding up the gelation process^. This effect, while certainly present, 
is however far too small to justify by itself the rapid fall in the gelation times for increasing salt 
concentrations, as observed in Figs. and Coupled with the fact that experimental ^-potential 
measurements have been observed to be relatively insensitive to counterion adsorption on the sur¬ 
face!^, we proceeded to approximate the surface potential to be constant, and set it equal to its 
dilute value of -30 mV^^. The second approach was therefore to proceed with the surface potential 
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Figure 5: The hydration force constant Fq diminishes with increasing salt concentration. These 
values are obtained from samples with a volume fraction (^udox of 0.133. 


fixed, so that we could fit the experimental data in Figs.andto obtain the hydration force con¬ 
stants Fq for each salt, and for each concentration. This is shown in Figs. and as function of 
salt concentration and particle volume fraction. It is clear from this framework that the reduction 
in Fq, and by extension in the repulsive hydration potential, upon increasing the salt concentration, 
is what controls the power-law relationship observed in Figs. and 



Figure 6: The hydration force constant Fq has a weak but systematic dependence on the particle 
volume fraction (^udox- 


The hydration force magnitude values Fq, are obtained at each salt concentration (c) from a 
one-parameter fit to the experimental data points in Fig. The resulting functions Fq{^{c)) are 
shown in Fig. (five salts) and Fig. (four different values of the particle volume fraction <|>Ludox)- 
Three trends are visible: (1) Fo(c) is higher for KCl than for any of the other ions, this means a non- 
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monotonic behaviour in terms of ionic size; (2) Fo(c) increases as Ludox concentration decreases; 
(3) Fo(c) decreases as salt concentration increases. If our effective potential is valid, we have to 
assume all three effects are related to the association of the counterions with silica surface charge 
groups (or else, other factors might be contributing to the interaction, and are being assimilated 
into these Fq trends). 

A plausible expected behaviour is that it is more difficult for the smaller, more hydrated mono¬ 
valent cations to approach and thus associate with the surface hydroxyl groups. Counterions ad¬ 
sorbed on the silica surface act as a repulsive force between particles. This explains why the 
hydration potential V/j increases monotonically as a function of the hydration diameter (see Fig.|^. 
Note that since Vf, FqA^ (eq. ??), the trend of Fq (Figs. and [^, which is observation (1) above, 
is more complex. A greater particle volume fraction corresponds to a greater total surface area, 
which is consistent with diminishing counterion association per unit surface area, and slight reduc¬ 
tion of Fq, which is observation (2) above. Also consistent with this picture is the fact that as we 
increase the salt concentration, the proportion of surface charge groups remaining free for coun¬ 
terion association dwindles and the ion-specific values of Fq begin to converge. Fig. shows that 
the hydration potential increases as the hydration diameter A lengthens, with the ordering of the 
salts by the relative strength of repulsive hydration forces remaining the same over all distances. 
Therefore the greater the hydration diameter of the counterion, the more long-ranged the hydration 
force becomes, and the greater the potential barrier to gelation. As a key result, the gelation time 
increases together with the hydration diameter of the ions in the Hofmeister series. 

In Fig. the hydration potential is observed to be a short-ranged monotonically decreasing 
function of the salt concentration (observation (3) considered above). The precise origin of this 
effect is not obvious. In our simple picture of counterion adsorption we expect more ions to 
adsorb, the higher the bulk concentration. A possibility, similarly to what proposed inl^, is that one 
needs to consider a loss of hydration shell when many ions adsorb, with a corresponding decline 
in repulsion. A final point we should remind the reader is that the surface charge density in the 
DLVO terms has been kept constant; this is unlikely to be strictly correct, but it is very difficult 
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Figure 7: The salt speeifie repulsion is very important at short range, as shown in (a) by eomparing 
the total interpartiele potential (Vtot, solid line) to the potential (VdvlOj dashed line) without the 
hydration potential V^, and in (b) by plotting the hydration eontribution by itself. The potentials 
deerease with salt speeies in the order LiCl-NaCl-KCl-RbCl-CsCl. Salt eoneentration is equal to 
0.3M for both plots. 


to do otherwise with the data at hand^^. Also, there are no experimental teehniques to aeeurately 
evaluate the surfaee potential, and standard zeta-potential measurements are not adequate for this 
task. So the question of whether this partieular trend originates from some physieal force or 
change in conditions at the gap, or an external force such as the bulk osmotic pressure, remains to 
be properly addressed in future studies. 

Of particular interest to this work, Trompette and co-workers have described in a series of pa- 
persl^^^^the effect of NHJ and Na'*' on the stability of colloidal silica, and found that samples with 
NH 4 aggregated much faster than those with Na"*'. They ascribed this to the different degrees of 
hydration of the ions. For the experiments described in this paper, only LiCl, NaCl, KCl, RbCl and 
CsCl were studied since the ammonium salt ion is somewhat acidic, and the pH itself also affects 
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salt concentration (M) 


Figure 8: With increasing salt concentration, the hydration potential decreases. The interparticle 
separation is set here at the hydration diameter for each salt species. Solid lines are pVtoi and 
dashed lines are /3l4, same color scheme as in|^ 

stability of colloidal silica. 

Experimental papers investigating short-ranged hydration forces have classically employed 
Surface Force Apparatus (SEA), and more recently also Atomic Force Microscopy (AFM). These 
techniques can be applied to measure the forces between smooth solid surfaces, lipid bilayers and 
biomembrane surfaces^. Some of the most recent studies include the use of AFM to estimate 
short-ranged hydration forces induced by multivalent salts^*— i the investigation of ion-specific 
Hofmeister effects between planar single-crystal sapphird^ and the measurement of charge inver¬ 
sion as a function of pH and salt concentration. These studies are in agreementl^^FQ] ^ trends 

presented here, whilst in other conditions the Hofmeister series is in reverse ordei®^, but we have 
not found experiments that can be directly compared to our results. We also note that our modeling 
shows that variations in the surface potential are insignificant compared to the dominant short- 
ranged hydration forces at close surface to surface separations. Ion-specific double layer pressure 
has been calculated using the full nonlinear Poisson-Boltzmann equation with the addition of the 
ionic dispersion energy between the ions and the two interfaces!^ By treating the electrodynamic 
ionic dispersion potentials on the same non-linear level as electrostatic potential, the ion-specific 
Hofmeister effects are recovered. We believe that there is great scope to combine studies such as 
this with experimental approaches such as in this manuscript, to help isolate and calibrate the vari- 
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ous possible mechanisms underpinning the hydration forces and bridge the gap between molecular 
and macroscopic observation. 

In conclusion, systematic experiments were carried out quantifying the gelation kinetics as a 
function of monovalent salt species. By using a simple model for the kinetics of cluster aggregation 
to fit experimental data of gelation times, at different salt conditions, for the first time it was 
possible to extract the magnitude of a non-DVLO hydration repulsion, that has a range set by 
the solvated diameter of the counterion. The very simple non-DVLO term used here is obviously 
coarse graining the detailed molecular mechanisms (the ordered “rigid” water layers repelling each 
other, and the energy required to “squeeze” these away as the particles approach contact) and is an 
effective semi-empirical term. This approach is powerful because experiments can then be readily 
fitted by a single parameter, the amplitude of repulsion force (Fq, function of salt concentration). 
The latter decreases with increasing salt concentration due to de-solvation of the counterions upon 
adsorption on the surface. The key finding is that the hydration repulsion correlates positively 
with the chaotropic nature and the size of the cation species. This framework and the molecular- 
level mechanism proposed here can be used in the future to devise tunable gelation protocols of 
nanoparticles by choosing the salt type. 

Methods 

All samples studied consist of colloidal silica, water and salt in different concentrations. Commer¬ 
cial silica colloids of Ludox HS-30, (Sigma-Aldrich), were used. Ludox HS size measurements 
vary from 16.7 nm^^to 18.5 nm® diameter, and we have taken 17 nm diameter as a value in our 
calculations. 

The Ludox was filtered prior to use, using Millipore Millex GS 0.22 jim filters, to ensure 
the removal of larger aggregates. The resulting Ludox “stock” used in most of the experiments 
had a density of 1.23 g/ml and contained 31.5% silica by weight. The pW of the original Ludox 
suspension (before filtration and dilution) was 9.8 (determined by the manufacturer), and was not 
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regulated during the experiments. Stoek solutions of salt were prepared using LiCl (99.0%, Sigma 
Aldrieh and 99%, Aeros), NaCl (99.0%, Sigma Aldrieh and analysis grade, Merek), KCl (99.9%, 
Fiseher and 99+%, Aeros), RbCl (99.8+%, Aeros) and CsCl (99+%, Aeros). All water used for 
preparation of stoek solutions and samples was of Millipore grade. 

The salt solution was always added last to make the final sample, of 1 ml volume. All samples 
were mixed on a vortex mixer for around 10 seeonds immediately after the addition of salt stoek. 
The gelation time was determined on a maeroseopie seale, by gently inverting the vials to observe 
the presenee (or not) of flow. When there was no diseemable flow for around 1 seeond, the sample 
was judged to be gelled. This eriterion is somewhat arbitrary, but was strietly adhered to, so that 
results of different samples are eonsistent. Fig. [^e) shows photographs of a sample before and 
after its gelation time. 

Preliminary runs were earned out, so that the gelation time was approximately known. This 
ensured that subsequent samples were not inverted unneeessarily (we did not notiee in any ease 
eorrelations between the gelation time and the frequeney of inspeetion). Identieal samples pre¬ 
pared on different days did sometimes show different gelation times, whieh eould be due to small 
temperature ehanges or minor differenees in sample eomposition or preparation. These differenees 
were never very large (around 10%) and do not affeet the observed trends. 

Supporting Information Available: Details the eomplete expressions for the attraetive van 
der Waals potential and eleetrostatie interaetion terms used in this work. This material is available 
free of eharge via the Internet http://pubs.aes.org, 
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